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1. Introduction 

High-energy collisions of heavy nuclei result in the emission of up to several thousands of 
particles. Such large numbers reflect the complexity of the collisions; yet, they allow for a wealth 
of different tools to characterize the processes that take place, in particular to assess the properties 
of the medium that is created. Among the possible approaches, the presence of these many parti- 
cles calls in a natural way for statistical descriptions, as well as for observables involving several 
particles. Various studies are thus devoted to extracting correlations between pairs, triplets, and 
more generally ?i-tuples of particles. 

The purposes of these investigations are manifold, as are the underlying physics pictures: To 
mention only correlations in momentum space, there are studies focusing on identical particles that 
are close in momentum (and in position) space, searching for evidence of (anti)symmetrization of 
the corresponding wave-function Other, closely related studies deal with pairs of non-identical 
particles with close velocities, to estimate the difference between their emission times [^. An- 
other phenomenon investigated through studies of the correlations between emitted particles is the 
anisotropy in the transverse plane of the emission pattern induced by the finite impact parameter 
between the nuclei ("anisotropic flow") [||]. The extended transverse-momentum range that be- 
came available with collisions at RHIC allowed novel investigations, involving particles with high 
transverse momenta pj, with a view to studying jets. A first kind of such study is that of correla- 
tions in azimuth between pairs of particles with both high transverse momenta [Q], to look for the 
presence of structures, close or away in azimuth to a high-/7r reference ("trigger") particle, that 
resemble the jets seen in pp collisions. The next step consists of studies of correlations between 
a high- pt trigger and low-joj- "associated" particles [^, so as to characterize the response of the 
created medium to the propagation of a high-joj- parton (or, less plausibly, hadron). 

Whatever the specific aim of a correlation study, it boils down to a simple principle, namely 
to try to identify the trace of some genuine dynamical effect in the joint two-, three-, M-particle 
distributions. To accomplish this, one needs to determine properly what the expectation would 
be for these distributions in the absence of such a dynamical effect. Now, such a "no-dynamics" 
M-particle distribution is not merely the product of M single-particle distributions, for there exists 
a trivial correlation between arbitrary final-state particles, due to the conservation of total mo- 
mentum, which imposes some constraints on the joint distributions [^]. The computation of these 
constraints will be discussed in section ^, using the general formalism from probability theory 
introduced in section ^ In section ^, I shall further discuss the meaning of this ever-present cor- 
relation due to global momentum conservation in the collision, and speculate on possible ways to 
take it into account in correlation studies. 

2. Probability distributions and cumulants 

Consider a collision with a total of N particles in the final state' (throughout this paper, A'^ is 
assumed to be large). The basic observable in studies of the momentum correlations between M 
particles among the N is the joint distribution dpj • • • dp^- To avoid normalization issues, it 

includes all particles, including the non-detected ones. 
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is more convenient to consider the joint M-particle probability distribution /(pi , . . . ,Pm)> which is 
by definition normalized to unity, and therefore (roughly) independent of the system size. 

By definition, M particles with momenta pi,...,pM are statistically independent from each 
other if and only if the corresponding joint probability distribution can be factorized into the product 
of the M single-particle probability distributions: /(pi , . . . ,pm) = /(pi) • • '/(Pm)- Reciprocally, 
if they are not independent, this factorization no longer holds and the joint probability distribution 
involves further terms: the joint probability distribution can be expanded into a sum over all prod- 
ucts of cumulants corresponding to distinct partitions of the M particles. For instance, at the two- 
and three-particle levels: 

/(Pl,P2) = /c(Pl)/c(P2)+/r(Pl,P2), (2.1) 
/(Pl,P2,P3) = /c(Pl)/c(P2)/c(P3)+/c(Pl)/c(P2,P3)+/c(P2)/c(Pl,P3)+/c(P3)/c(Pl,P2) 

+ /c(Pl,P2,P3), (2.2) 



where the single-particle cumulant is equal to the single-particle probability function, /c(p) = 
/(p). The M-particle cumulant /c(pi, . . . ,PAf) corresponds to the "genuine" correlation between 
the M particles. The physical interpretation of the cumulant expansion is straightforward: the 
joint M-particle distribution depends not only on the genuine M-particle correlation, but also on 
all the possible correlations involving subsets among the M particles. As an example, think of 
two particles emitted exactly back-to-back with p2 = — Pi (as the two pions from a decaying p 
meson in the rest frame of the latter). The two-particle probability distribution reads /(pi,P2) = 
/(pi)5(pi +P2) (where, in the p nn case, the single-particle distribution /(pi) actually reduces 
to the angular distribution since |pi| is also fixed), while the corresponding cumulant is obviously 
non- vanishing, see equation (2.1). Other illustrations of the difference between the distributions / 
and the cumulants fc can be found in reference [^. 

While equations {2A)-('L2), and so on, can be inverted one after the other to yield the cumu- 
lants as functions of the joint probability distributions, there is a more systematic way to perform 
the same operation. First, one defines a generating function of the joint multiparticle probability 
distributions: 



G{xi,...,xn) = l+;ci/(pi)+X2/(p2)H ^xiX2f{puP2) + ---, (2.3) 

and similarly for every order M, where xi,...,xn are auxiliary variables. Given this generating 
function of the joint probability distributions, the function that generates the cumulants is simply 
its logarithm [||] : 

lnG(;ci, . ..,xn)= xi/c(pi) +X2/c(p2) H hxiX2/c(pi ,P2) + • • • . (2.4) 

Thus, the knowledge of one of these functions automatically translates into that of the other. 
In the following section, I shall use this property to sketch the computation of the multipar- 
ticle cumulants due to the momentum-conservation constraint, starting from the expression of 
the joint probability distribution for M particles. Moreover, I shall also use "scaled" cumulants 
(pi, . . . ,Pm) = /c(pi, • • • ,Pm)/[/(Pi) • • -/(Pm)]. 



3 



Multiparticle correlations and momentum conservation 



Nicolas Borghini 



3. Multiparticle cumulants from momentum conservation 

As stated in the introduction, the purpose is to determine how global momentum conservation 
affects the joint multiparticle probability distributions of final-state particles in a large-multiplicity 
event like a heavy-ion collision. That is, given A'^ particles with momenta pi, . . . , pat obeying 

the constraint pi H hpv = 0, what is the resulting M-particle cumulant? This question was to 

my knowledge first addressed (albeit semi-quantitatively) in the two-particle case in reference [^]. 
A quantitative estimate of the two-particle correlation was then derived in the frame of some 
anisotropic-flow measurement ||^, then independently rediscovered in the same context [|T^. In 
both these cases, the computation relied on the use of the central-limit theorem for the distribution 
of the sum of « 1 (yet n < N) uncorrelated momenta. Eventually, a general approach to com- 
pute the cumulants to arbitrary order was introduced in reference [[III], making use of the generating 
functions of joint probabiUty distributions and of cumulants and of a saddle-point integration. 

The starting point of these calculations is the expression of the joint M-particle probability 
distribution under the momentum-conservation constraint: 

/ M \ „ N 

n^(P7) /5(pi+---+Piv) n [^(Pi)dP;] 

/(Pi,---,Pm) = jj , (3.1) 

/5(pi + --- + pw)n[F(p,-)dpy] 

where F(p) is the single-particle probability distribution "unrenormalized" for the momentum con- 
servation constraint: to leading order in 1 /N, it equals the measured single-particle probability 



distribution /(p), yet they actually differ at the next-to-leading order []1Q|]. In the previous calcula- 
tions [^-|l|], F was assumed to be isotropic; here we shall relax this assumption and consider the 
more realistic case of collisions with anisotropic expansion (flow) in the transverse plane.^ Quite 
obviously, the Dirac 5 in equation ( p?T| ) represents the constraint from global momentum conserva- 
tion: in its absence, /^(p) = /(p) and the joint probability distribution would simply factorize into 
the product f (pi) • ■■F{t^m)- 

To derive the cumulants arising from global momentum conservation, the most convenient and 
systematic way is to introduce equation ( p| ) in the expression of the generating function ( p^ ) so 
as to compute the latter [pl|]. After introducing a Fourier representation of the Dirac distribution, 
one finds that G{x\ ,xn) can be expressed as the integral over the Fourier conjugate variable k of 
the exponential of a function NJ^{k), which also depends on the auxiliary variables xj.^ Since N is 
supposed to be large, this integral can be performed by a saddle-point approximation, provided one 
finds the position ko of the maximum. The key to obtaining the successive cumulants is then first 
to solve to a given order in x/N the equation giving ko (a solution to the order M — 1 is required for 



^The recipe for extending tlie calculation to the non-isotropic case was briefly given in reference [jll]], and the 
corresponding expression for two-particle correlations can be found in reference 

■^More precisely, ^ only depends on xj through combinations XjF{pj)/N: this allows one to replace Xj in G by 
Xj = XjF(p j) — which amounts, to leading order in 1 /N, to replacing the cumulants fc by the scaled cumulants fc — 
and to derive the scaling with A' of the cumulants nil| : the M-particle cumulant scales as 1 /N^^ ' . 
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the M-particle cumulant): 



where the angular brackets denotes a F(p)-weighted average. Then, one uses this expression of ko 
to compute the value of A'^^(ko): 

A^^(ko) -A^ (ln(e^''°-^) + I • (3-3) 

The coefficient of Xj^ . . .xj^ in A'^^(ko) is then the scaled cumulant /c(pi , . . . ,Pm)- For instance, 
one finds (we assume for the sake of brevity that (p) = 0, otherwise one only need replace p by 
p' = p — (p) in the following formulas) 

f Cn n ^ - P^^^P^,^ PhyP2,y Pl,zP2,z 

/c(Pl,P2j — — ^r/.-7\ T77ZT\ ^TTZtT' '3-4) 

/c(Pl,P2,P3) = - 




(3.5) 

In these expressions, the x-, y- and z-directions are the principal axes that diagonalize the tensor 
(pCgip): in practice, one axis will lie along the beam direction, one along the impact parameter of 
the nucleus-nucleus collision, which is the only preferred direction in the transverse plane, and the 
third one will be perpendicular to the other two. If F{p) is isotropic, {p"^) = {pi) = {pi), so that 



one recovers the formulas given in reference 1 11]. 



In nucleus-nucleus collisions at ultra-relativistic energies, the mean square momentum along 
the beam direction (z) is typically significantly larger than the mean square transverse components. 
As a consequence, the terms involving the z-component in the expressions of the cumulants are 
much smaller than those involving the other two components, especially when one considers final- 
state particles emitted close to mid-rapidity. Therefore, I shall drop these longitudinal-momentum 
terms from now on, which amounts to considering only the constraint imposed by the conservation 
of total transverse momentum pj-. 

4. Defining a minimally-biased background for correlation studies 

The behaviour of the three-particle cumulant resulting from the momentum-conservation con- 



straint, equation (pj), was discussed in some detail in reference [13| in the case of an isotropic 



transverse emission of particles. Here, I shall now discuss further the expression of the two-particle 
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cumulant and its content. The physical meaning of equation (3.4) is intuitive: momentum conser- 
vation induces an "anti-correlation" between the momenta of final-state particles. This is of course 
reminiscent of the situation in which only two particles are emitted, with back-to-back equal mo- 
menta (note, however, that the calculation reported here cannot cover that case, since it assumes 
A'^ » 1). Additionally, the anti-correlation is largest between particles with larger momenta; it is 
also largest in events with a smaller total multiplicity A'^. 

Let me emphasize the implications of the expression of the two-particle cumulant. In colli- 



sions in which particles are in the final state, equation ( [3.4| ) implies that the joint two-particle 
probability distribution reads (assuming first that particles are emitted isotropically, which is a fair 
approximation in the transverse plane for central collisions, so that {p\) = {p^) = {pj) /2): 

/(Pn,Pr2) =/(Pn)/(Pr2) (^-^T^^^) • ^^.l) 

In other words, although the particle emission is isotropic, i.e. /(pr) depends on pj only, yet 
given a "trigger" particle with transverse momentum pj-j, the conditional probability to find an 
"associated" particle with transverse momentum ipj2 is not isotropic: 

/(Pri,Pr2) ^A 2p7-i-pr2 



/(P„|Pn) = =/(p.,.) ^. -i^^j #/(Pn). («, 

Thus, the conservation of transverse momentum induces a sinusoidal modulation of the probability 
distribution of "associated" particles, with a minimum of the distribution along the trigger-particle 
momentum: there is a larger probability that the momentum of the associated particle points in the 
hemisphere opposite to that of the trigger particle. The amplitude of the modulation increases with 
the values of both the trigger momentum pn and the associated momentum pt2, and it decreases 
with increasing A'^. 



One recognizes in equation (4.2) the difference between marginal [/(pr2)] ^^i^ conditional 
[/(p7'2|p7'i)] probability distributions [|8]]. Since the former is — up to a normalization factor — 
the measured single-particle distribution, it is most tempting to use it in studies of two-particle 
correlations; yet one must rather use the latter once a first particle momentum pj^ in the event is 
fixed as a reference. This is not to be unexpected: by removing from the event this trigger particle, 
one obtains a collection of final-state particles (including the non-measured ones) which is not a 
"valid" event, since the sum of the transverse momenta of the particles does not vanish. Hence there 
is no reason why the single-particle probability distribution for this collection should be the same 
as for real events satisfying momentum conservation."* Thus, in studies of two-particle correlations 
in which the transverse momentum p^j of one of the particles has been fixed — thereby implic- 
itly selecting a subset of the whole available sample of events — one should consider the condi- 



tional probability distribution /(pr2 Pri) of equation (^.2), including the momentum-conservation 



constraint, as the proper "reference" distribution of associated-particle momenta pr2. over which 



'^In fact, from the mathematical point of view, these collections of A' — 1 particles with non- vanishing total transverse 
momentum J^Pt" = — pj- [ are equivalent to events with a mean first-harmonic anisotropic-flow component (directed flow) 
vi = —pTilYJ^^ 1)(P7')]> the direction of the equivalent reaction plane being that of pj-j . Therefore, the corresponding 
single-particle distribution includes a 2v\ cos((p — (p\) azimuthal modulation, which is absent from the original single- 
particle distribution of the real events, which we assumed isotropic. 
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dynamical effects are to be investigated. This is admittedly slightly unsatisfactory, since the con- 
ditional probability distribution is not measured, and depends on two unknown quantities: the 
total number N of final-state particles and the mean square transverse momentum (pj-) (albeit only 
through their product).^ Yet they can be estimated, as was done to take into account the conserva- 



tion of total momentum in analyses of anisotropic flow [ |10| , |14| ] or in femtoscopy studies [|12|]. Only 
at this price can one determine the "minimally-biased background" which is what the distribution 
of associated particles would look like in the absence of non-trivial correlations. For Au-Au colli- 
sions at RHIC energies, the corresponding correction might be of at most 1 or 2% when choosing 
a high-pj trigger particle; yet this is of the same relative magnitude as the effects that are mea- 
sured iQ], so that such a precision is necessary if one wants to establish the existence of specific 
dynamical effects and to quantify their importance. 

Let me now comment on the two-particle cumulants in the more general case where particles 
aie not emitted isotropically in the plane transverse to the beam. More precisely, I shall consider 
particles emitted with a mean second-harmonic transverse anisotropy (elliptic flow) V2 defined by 



V2 = {pI — Py)/{pl + Py).^ This definition yields at once the identities {pD = (1 +V2)(;?f-)/2 and 



(Pv) = (1 ~'^2){Pt)/'^' which one can insert in equation (p^: 



from which one deduces the conditional probability distribution /(Pt2|p7'i) of particles associated 
to a trigger particle with transverse momentum pn- One sees that the effect of the constraint 
from momentum conservation is larger in the j-direction, i.e. perpendicular to the nucleus-nucleus 
impact parameter, than in the x-direction. That is quite normal, since particles with a transverse 
momentum along y are less numerous than those pointing in the x-direction (remember that values 
of V2{pt) =0.16 were reported for pt>2 GeV/c in minimum-bias Au-Au collisions at RHIC [p^, 
meaning that twice more charged hadrons are emitted in the x-direction than perpendicular to it). 
As a consequence, if the trigger particle is chosen with pj i in the j-direction, only few particles 
are present to balance its momentum, so that those few ones are more strongly correlated to it than 
if we were considering a trigger particle with a momentum in the x-direction. 



The dependence of the strength of the two-particle cumulant ( |4.3[ ) on the azimuths of the two 
particles means that the correction for the momentum-conservation effect has to be performed with 
some care. Let me illustrate that with an example. I have already argued above why one should use 
the conditional probability distribution /(pr2|Pri) instead of the marginal one when investigating 
the possible structures associated with the presence of a (high-pj) trigger particle. If only the 
average correction, corresponding to the isotropic case (p~T|), were used instead of the azimuthally- 
dependent one (p~3|), this would yield an over-correction (resp. an under-correction) of /(pr2|pri) 
for trigger particles emitted along the nucleus-nucleus impact parameter (resp. emitted along the y- 
axis). As a result, one would observe a small (using typical RHIC values for A^, (pj-) and V2, of order 



^Furthermore, the averages (• • •) involve the non-measurable distribution "in the absence of momentum conserva- 
tion" F(p), rather than the physical distribution /(p) Jl^. This is however but a minor issue, as the difference between 
F and / is of subleading order in 1 /N, while here we only consider leading-order quantities. 

^This definition yields values of V2 that are typically a factor 2 larger than those obtained with the more conventional 
definition V2 = ((pJ ^ Pv)/(.Px + Pv))- s° '^^at at RHIC energies V2 ~ 0.1 in mid-central Au-Au collisions. 
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< 0.005) spurious bump (resp. dip) at 180" away from the trigger in the conditional probability 
distribution for t^j2- This azimuthally-dependent spurious structure due to an inaccurate definition 
of the "minimally-biased background", which is to be the reference over which correlations of 
dynamical origin can be observed, would prevent any accurate determination of these interesting 
correlations, by mimicking irrelevant features. Further implications of the azimuthal dependence 
of the two-particle cumulant due to momentum conservation are discussed in reference [[l6|]. 



In summary, I have recalled that the general purpose of studying correlations is to yield ev- 
idence of phenomena that go beyond trivial expectations. In the specific context of high-energy 
collisions, correlations between any number of final-state particles are induced by the conservation 
of total momentum, which are not of dynamical origin. These uninteresting correlations can be 
computed, thereby allowing one to define a "minimally-biased background", including the effect 
of total-momentum conservation, which is viewed as the reference over which genuine dynamical 
effects might be revealed. 
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